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We study the dynamical equilibrium correlation function of the polaron- 
dressed tunneling operator in the dissipative two-state system. Unlike the po- 
sition operator, this coherence operator acts in the full system-plus-reservoir 
space. We calculate the relevant modified influence functional and present the 
exact formal expression for the coherence correlations in the form of a series 
in the number of tunneling events. For an Ohmic spectral density with the 
particular damping strength K = ^, the series is summed in analytic form 
for all times and for arbitrary values of temperature and bias. Using a dia- 
grammatic approach, we find the long-time dynamics in the regime K < 1. 
In general, the coherence correlations decay algebraically as t~ 2K at T = 0. 
This implies that the linear static susceptibility diverges for K < ^ as T — > 0, 
whereas it stays finite for K > ^ in this limit. The qualitative differences with 
respect to the asymptotic behavior of the position correlations are explained. 



I. INTRODUCTION 

The simplest model that allows to study the interplay of tunneling and dissipation is the 
spin-boson model . Despite its simplicity, it exhibits generic features of many complex 
systems and has found widespread applications in physics and chemistry. It has been adopted 
to describe diverse systems, like the tunneling of atoms between a surface and the tip of an 
atomic- force microscope ||, or the dynamics of the trapped flux in a quantum interference 
device [|J , to mention a few. 

For Ohmic dissipation, the spin-boson model shows a transition between coherent and 
incoherent tunneling where the location of the transition depends on the damping strength 
and on the bias 0. Most investigations have been done for the non-equilibrium expecta- 
tion value (a z (t)), representing the population difference between the two localized states, 
and for the a z autocorrelation function, describing position or population correlations ||[| . 
For the latter quantity, the analogy with the Kondo model and the 1/r 2 Ising model has 
been utilized in numerical computations. For Ohmic damping and zero temperature 
and bias, (a z (t)) shows a transition between damped oscillations and incoherent relaxation 
exactly at K = \ [|l(J . Recently, it has been argued fll] that the quality factor of the oscil- 
lation is exactly given by Q = cot[7ri ; r/2(l — K)], which again gives K = | for the transition. 
For the antisymmetrized a z autocorrelation function, the same quality factor and thus the 
same transition point was found numerically [[| . This is not in contrast to the earlier result 
K = ^ |8|JT^1, since there a different criterion for the transition was applied ||[13 . 



It has been shown that the particular initial preparation plays a crucial role for the 
long-time behavior at zero temperature. The factorized system-reservoir initial state for the 
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expectation value (& z (t)) leads to exponential decay flO| . |TT|j , whereas the symmetrized a 



equilibrium correlation function with a correlated initial state decays algebraically as 1/t 2 
for K < 1, as shown for the spin-boson model || and the related 1/r 2 Ising and fermionic 
models • The power 2 in the algebraic decay law is a signature of Ohmic dissipation. 

Recently, focus has been put on expectation values [rjj] and equilibrium correlations fl"5| 
connected with the tunneling operator a x . The expectation value (a x (t)) and the equilibrium 
autocorrelation function of a x have been found to be nonuniversal, i.e., vanish in the scaling 
limit |14],16|1. Here, we study the equilibrium autocorrelation function of a polaron-dressed 



tunneling operator a x which includes the adiabatic dynamics of the bath modes P],|T7[ . This 
function is universal and measures correlations of the off-diagonal elements (coherences) of 
the density matrix. We present the exact formal solution for the coherence correlations in the 
form of a series in the number of tunneling events. We then analyze the resulting expression 
in various limits. In particular, we work out the differences in the asymptotic decay between 
the position correlations and the coherence correlations. Our analytical real-time approach 
is complementary to the recent imaginary-time numerical studies in Ref. ||15|| . 

In Section II, we formulate the problem and introduce the correlation function C x (t) of 
the coherence operator a x . Since this operator acts both in the system and bath space, the 
elimination of the bath modes has to be reconsidered. The relevant considerations leading 
to a modified influence functional are given in Section III. These results are used in Section 
IV to determine exact formal expressions for C x (t). In Section V, we present the analytical 
solution for C x (t) for the special value K = ~. Section VI is devoted to the regime K = | — n 
with k C 1. Finally, in Section VII we show that the asymptotic decay of C x (t) at zero 
temperature is algebraic with a i^-dependent power for < K < 1. 



II. FORMULATION OF THE PROBLEM 

It has been well established that the dissipative dynamics of a particle in a double-well 
potential can effectively be described at very low T by the spin-boson model 
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Here the basis is formed by the two localized states \R) and \L) of the double- well which 
are eigenstates of a z with eigenvalues +1 and —1, respectively. The position operator is 
q = ao z /2 with a z = \R)(R\ - \L)(L\. The tunneling operator a x = \R)(L\ + \L)(R\ 
transfers the particle between the two wells with tunneling amplitude A. The second term 
in Eq. ( |2.2| ) describes an externally applied bias energy e. The effect of the thermal bath on 
the system's dynamics is included in the spectral density 
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The important case of an Ohmic bath is described by 
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where r] is the viscosity, K is the appropriate dimensionless damping strength, and oo c is a 
cut-off for the bath modes. We are interested in the regime A <C u c , in which A and uo 
form a renormalized frequency scale ffl 



c 



A r = A(A/u c ) K ^- K ^ . (2.5) 

A quantity is called universal, if it is a function of A r alone, i.e., there is no other uo c 
dependence than given by Eq. (|2.5|) . Vice- versa, any extra dependence on u c is non-universal: 
performing the scaling limit A r /uj c — > with A r fixed, this contribution vanishes. Both 
the mean value of a z and its equilibrium autocorrelation function are universal [l],f|- On 
the other hand, the expectation value (a x (t)) is equipped with an overall factor A r /A = 
(A r /u} c ) K and therefore vanishes in the scaling limit |TJ| . 



Here we concentrate on the a x equilibrium autocorrelation function. As observed in 
Ref. ]n| and explained in the sequel, the equilibrium correlation function of the bare a x does 



not satisfy the above universality criterion. To overcome this shortcoming, we consider a 
modified tunneling operator which takes into account the adiabatic displacement of the bath 
modes during the tunneling process. The transformation to a basis of displaced harmonic 
oscillators states is accomplished by the polaron unitary transformation M 



U = exp {— ia z Q/2h} 
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The set of displacements is given by {s a } = {ac a /m a uj'^}. The polaron transformed tunnel- 
ing operator a x = Ua x U~ l reads 

a x = \R) (L\ exp{-in/h} + \L) (R\ exp{in/h} (2.7) 
= \R)(L\ / dx \x) (x-s\ + h. c. , (2.8) 



where we have introduced the compact notation 

s = {s a } , x = {x a } , J dx = Y\ J dx a . (2.9) 

a 

The bare a x acts in the Hilbert space of the two-state system (TSS) alone, whereas the 
dressed operator a x acts in the full system-plus-reservoir space. From the coordinate repre- 
sentation ( |2.8| ), we immediately see that the operation of a x transfers the particle from one 
localized state to the other and simultaneously shifts each bath-oscillator by the displace- 
ment ±s a ("polaronic cloud"). In this sense, a x generates coherent tunneling between the 



two localized states and can be called coherence operator [18]. The coherence correlations 



are then described by the equilibrium correlation function of a x , 

C x {t) = {a x (t)a x (0)) p = Tr [^(^(0)^] , (2.10) 

where Wp = e - ^/ Tr [e _/3H ] is the equilibrium density matrix of the global system, and 
a x (t) is the Heisenberg-representation of a x with respect to the untransformed Hamiltonian 
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( j2.ip (cf. Ref. [0). The associated response function Xx(t) = ( — Q(t) lmC x (t) describes 
the linear response of the system to a coherence inducing perturbation H peTt oc a x . 

It is convenient to consider C x (t) as the mean value of & x (t) with respect to the "density 
matrix" W = a x (0)W@. Switching to the Schrodinger-picture, Eq. (|2.10| ) becomes 

C x (t) = r Fr[a x W(t)], (2.11) 

where the time- dependent "density matrix" W(t) = exp(—iHt/h)W(0) exp(iHt/h) obeys 
the initial condition W(Q) = a x Wp. Inserting the expression fl2.8Q for a x into Eq. (p. 11] ) and 
performing the trace, we find that C x (t) is the sum of the off-diagonal matrix elements 

C x (t) = p[ s ^(t) + p { - s l(t) ■ (2.12) 

Since a x acts also in the bath space, p[ s J (t) is different from the usual reduced density matrix 

as it appears, e.g., in the a z correlation function 0. We shall refer to pf}{t) as the "shifted 
reduced density matrix" (SRDM). To be general, we now give the discussion for a continuous 
variable q, and return to the two-state system only in Section IV. We have 



P W (9/>9/>*) = J dx f (lf, x f + s \W(t)\q' f ,x f ) 

= J dqidq^dxf dxidx' i K(qf,Xf + s,t;qi,Xi,0) (2-13) 
x K*(q' f , x f ,t; q' i: x' i ,0)(q i , xi\ W(0) , 

where K(qf, Xf + s, t; q i} Xi, 0) is the usual Feynman propagator that may be expressed as a 
path- integral. The matrix elements of the initial "density matrix" W(0) = o x W$ read 

(q^XilW^lq^x't) = (-q^Xi-ssgniqdlWpl&x'i) , (2.14) 

where the bath coordinates are also affected by the preparation. Now, it remains to integrate 
out the shifted bath degrees of freedom in the expression (|2.13|) . 



III. GENERAL INITIAL PREPARATIONS 



The standard Feynman- Vernon approach which may be used to eliminate the bath de- 
grees of freedom relies on the assumption of a factorized system-bath initial state. For 
ergodic systems, it is possible to obtain a real-time description also for a thermal initial 
state 0. This approach can be generalized to special classes of correlated initial states by 
introducing a preparation function pOf . For the case of C x (t), the method in Ref. |2(J has to 
be reconsidered, since the initial preparation also involves the bath. To proceed, we define 
a generalized preparation function A G (gj,^; q, x,x') by 



(qi,Xi\ W(0) \q'i, x[) = Jdq dq' dx dx A G (gi,g-; q,q'; Xi,x[; x,x') (q,x\ Wp \q',x') . (3.1) 
Comparing this form with Eq. ( |2.14j ), we see that the preparation function factorizes as 



A G (?i,9i; q,q';%h%i, ^0 = A s (%,?-; ^gOVO*^; x , x ') 



(3.2) 
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where the system's and reservoir's preparation functions are given by 

A R (x«, x^; x, x') = 5(xj — x — s sgn(gi)) <5(x^ — x') . 
With the form ( [OP , the evolution of the SRDM ( [2~T3D is given by 



P {s \qf,q'f,t) = j dqidq'idqdq' J G (qf,q' f ,t;q i ,q' i ;q,q')X s (q i ,q' i ; q,q') , 
where the generalized propagating function reads 

■/g(?/j q'fi t'1 Qu q'i'i 9) q') — J dxfdxidx' i dxdx'K(qf,Xf + s,t;qi,Xi,0) 

x (q, x\ Wp \q, x ) K*(q f ,x f , t; q i} x iy 0) A R (xj, x-; x, x 



(3.3) 
(3.4) 



(3.5) 



(3.6) 



For an ergodic system, the thermal density matrix Wp in Eq. (|3.6|) can be expressed as 
follows. Describe the global system at a time to < by a factorized density matrix, the 
system being in a position eigenstate, say \qo), and the reservoir being in thermal equilibrium, 
W(t ) = \q )(q \ ® e~^ R /Tr [ e _/3HR ] @, and let it evolve out of this state under the full 
Hamiltonian. Then, if t is sent to the infinite past, the system has reached at time zero the 
correlated initial state (q,x\ Wp \q',x'). With these considerations, we may rewrite (3.6) as 

Ja(qf, q'f, *; ?*> e<> 0+ ; ^ °~; ?o, so> *o) 
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where iS s [<?] is the action corresponding to the system Hamiltonian 
integrations are over all paths q(t') and q'(t') which satisfy the constraints 



•J~ g[?) q'\ s ] > (3-7) 
. The functional 
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q(0 + ) 
q'(0+) 
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q(t) = q f (3.8) 
q\t) = q' f . (3.9) 



All the effects of the bath onto the system's dynamics are captured by the generalized 
influence functional 



(3.10) 
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where <S Ri i[x, q] is the action corresponding to the reservoir and interaction terms in Eq. Q2.1| ). 
The paths x(t') and x'(t') are subject to the constraints 
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x(t) = x f + s (3.11) 
x'(t) = x f . (3.12) 



Compared to the standard Feynman- Vernon influence functional, there are two differences. 
First, the endpoint of the x-path is shifted by the displacement s, i.e., the bath does not end 
up in a diagonal state at time t. Second, the reservoir paths x(t') and x'(t') are discontinous 
at time zero, depending on the reservoir's preparation function A R (xj, x\\ x, x'). 
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IV. EXACT FORMAL SOLUTION 



Having obtained an explicit expression for the SRDM at time t for a general initial 
preparation, we can now write down the exact formal solution for the coherence correlations 
C x (t), Eq. (|2.12|) . Inserting the preparation function ( |3.3|) into the SRDM (|3.5|) , we get 



C x (t) = lim VI J g(<?/, ~Qf, t; q h g-,0 + ; -q h g-,0 ; g , q ,t ) 



t(j— *— 00 



where g i; g 4 ', g/ = ±a/2. Note that there is a jump in the g— path at time zero. Inserting 
the reservoir preparation function ( p.4|) into the influence functional ( |3.10|) , the x'(t')-path 



( |3.12| ) turns out continuous at time t' = 0, whereas the x(t')-path is discontinuous. 



Because of the integration over 2j, the constraints (|3.11|) may equivalently be expressed as 



x(t ) = x x(0 ) = Xi x(0 + ) = Xi + s sgn(gj) x(t) = Xf + s sgn(g/) . (4.1) 

Consider first the contributions to C x (t) with qi = qf. In this case, the shifts in the 
x-path at times t' = + and t' = t are equal, and thus we can eliminate the shift at positive 
times by defining modified reservoir coordinates according to 

x(t') =x(t')-ssgn( qi )Q(t') . (4.2) 

The path x(t') is continuous at t' = and obeys x(t) = Xf. In the shifted coordinate, the 
bath is in a diagonal state at time t. As the action 5" RI [x, q] appearing in the influence 
functional is quadratic both in x(t') and in q(t') and bilinear in the coupling, the second 
term in Eq. ( |4.2|) can be absorbed into a modified g-path which is continuous at t' = 0, 

q(t') = q(t')-2 qi ®(t') ■ (4-3) 

Writing the influence functional in terms of the paths q'(t'), x'(t') and the modified paths 
q(t'), x(t'), the displacement s is completely eliminated from the description. Thus, after 
integrating out the bath degrees of freedom, we end up with an influence functional that is 
of the standard Feynman- Vernon form for a factorized initial state at time to, 

^ G [q,q' ] ssgn(q t )]=^[q,q'] . (4.4) 

All effects in a x induced by the polaronic cloud are in the modified path q{t'). 

Next, consider the contributions to (^4.1| ) with qi = —qf. Now, it is not possible to end 
up with an influence functional of the form ( (4.4p in which the shifts of the bath modes are 
fully absorbed into a modified path q(t r ). In the usual charge picture (see below), the case 
qi = —qf corresponds to sequences of charges that violate overall neutrality. As a result, 
A and uj c cannot be combined to a function of A r alone. Instead, each contribution comes 
with an additional factor (Ar/iOc)^ and therefore is nonuniversal. Thus in the scaling limit 
A r /u c — > 0, all contributions with g^ = — qf vanish. 

With the above, the correlation function is now given by 

C x (t) = lim V J G (q h -q h t; q h q[, + ; -g i; g-, 0"; g , g , t ) . (4.5) 

to—*— oo ' 4 
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At this stage, it is important to note that the free propagators in the propagating function 
depend on the original paths q(t') and q'(t'). The concept of modified paths is only used to 
express the generalized influence functional in the standard Feynman- Vernon form. For the 
evaluation of Eq. (|4.5| ), it is convenient to introduce the linear combinations 



V(t)= [q(t) + q'(t)]/a 
m=[l(t)-q'(t)}/a 



(4.6) 



describing propagation along the diagonal of the density matrix and off-diagonal excursions, 
respectively. For the two-state system, these paths are piecewise constant with jumps at 
times tj. As usual, the time intervals t 2 j < if < ^j+i m which the system is in a diagonal 
state are called sojourns, while the time intervals £2^-1 <t'< t 2 j spent in an off-diagonal 
state are referred to as blips. A sojourn is labelled by rjj = ±1, depending on whether the 
system is in state RR or LL. Similarly, £j = ±1 describes a blip in which the system is in 
state RL or LR. The lengths of the sojourn- and blip-intervals are denoted by Sj = t 2 j + i — t 2 j 
and Tj = t 2 j 



£2,3-1 j respectively. All paths that contribute to the correlation function 



start out from the initial sojourn i] at time to an d end in the blip state £ at time t. According 
to their behavior at time zero, they can be divided into two groups. In group A, the system 
jumps at time zero from a sojourn to a blip state [q[ = —qi at time + ). In group B, the 
system hops at time zero from a blip to a sojourn state (g- = qi at time + ). A general path 
with n blips at negative and m blips at positive times can be parametrized by 



n+m— 1 

v(t')= £ %[©(*' -%)-e(t'-£ 2i+1 )] 

3=0 
n+m 

3=1 



(4.7) 



with £2n+2m = £ Plj ■ For group A, we have t2 n +i = 0, whereas t 2n = for group B. According 
to the boundary conditions, the paths are subject to the constraints 



(group A) 
(group B) 



Generic contributions to group A and group B are sketched in Fig. |T[ 



(4.8) 
(4.9) 
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t>2n+2m = t 
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n+l 



2/i 







£n+m| 

£2n+2m = t 



FIG. 1. Path contributions to group A (top) and group B (bottom). The steps represent blips 
of either sign, and sojourns are indicated by the baseline. 
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Thus the correlation function is built up by two parts that correspond to these two 
different path classes. We have C x {t) = C x {t) + C x {t) with 

Ctit) = lim V J G (£, t; C + ; r? = 0"; % , t ) , 

(4-10) 

C*(t) = lim Vj G (e,t;r ? = e,0 + ;-e,0-;r /o ,t ) . 

to— »— oo ' • 

The path sum is over all sequences of blips and sojourns and implies time-ordered integration 
over the jump times. We introduce the compact notation 

ft pt ftk + l + l ftk+3 pO ptl 

\ V kjl {tj} = / dt k+l+1 dt k+ i... I dt k+2 j dt k ... dti . (4.11) 

Jt Jo Jo Jo Jt Jt 

Here k and / represent the number of flips in the time regimes to < t' < and < t' < t, 
respectively. For group A, we have k = 2n, I = 2m — 2, whereas for group B k = 2n — 1 
and / = 2m — 1. Each transition in Eq. ( J4.ll ) comes with a factor ±zA/2. There are two 
additional transitions at time zero, t k+1 = 0, and at time t, t k+ i +2 = t. These two hops, 
however, come without a factor ±zA/2 since they are not dynamical. The jump at time zero 
is enforced by the operation of a x , whereas the jump at time t is introduced for convenience 
(cf. Ref. ||21|| ). The amplitude to stay in a sojourn is unity, while the amplitude to stay in 
blip £j is given by exp(ie^r ? ). Thus, a full path gives for both groups a factor 

-Tto^-A 2 /*)"-^- 1 ^,™ (4.12) 

with the bias term 

n+m 

D n>m = exp (it ^ £i r i) • ( 4 - 13 ) 
j=i 

Before discussing the modifications due to the polaron transformation, consider the stan- 
dard influence functional. Performing integrations by parts, it takes the form [^] 

T\m i\ = exp j jf dt'j^ dt" [i(t')S(t' - t")i(t") + ii(t')R(t' - t")r](t") ] } , (4.14) 

where the kernels S(t) and R(t) are the real and imaginary part of the second integral of 
the bath correlation function. In the limit u c t ^> 1, we have |2| 

S(t) = 2K\n [{Ufiujir) sinh(7rf/ft/3)] , (4.15) 
R(t) = irKsgn(t) . (4.16) 

Because of the form (|4.7| ), the velocities in Eq. ( 4.14| ) consist of a series of delta-functions 
centered at the flip times. This suggests to regard the blip- and sojourn-paths as sequences 
of charges: blip-charges interact with each other through the kernel S(t), while the sojourn- 
charges interact with the blip-charges via R(t). Substituting the paths ( |4.7| ) into Eq. Q4.14| ), 
the influence functional takes the form 

•F n,m G n rn H n m . (4.17) 



s 



The factor G n ^ m contains all the inter- and intrablip-interactions, 



G, 



n+m 



n+m j — 1 

^2 0& A i,fc 



A 



j,k 



«p | - Yl " E >Zvs/,-\ / ,,_ , (4.18) 

j=i i=2 k=i 

S2j,2k-1 + S , 2j-l,2fe — S-y^k — S2j'-l,2fc-l , (4.19) 

where S^g = S(t p — t q ). The sojourn-blip interactions are captured by the phase factor H n rn . 
With the form Q4.16Q , each sojourn only interacts with the subsequent blip, 



H, 



n,m 



exp 



n+m—l 
k=0 



k+1 



(4.20) 



Substituting Eq. (|4.15|) into the term A 2n+2m G n m , the quantities A and oj c are combined 
into a factor Ar 2X ^ n+m ) ; where A r is the renormalized tunneling frequency, Eq. ( |2.5| ). 
The autocorrelation function of the bare a x depends on the standard influence functional 
( |4.17| )- (|4.20| ). In this case, however, there appears the quantity A 2n+2m ~ 2 G n ^ m because 
the A-factors of the two blip-charges at time zero and time t are missing. Therefore the 
autocorrelation function of the bare a x comes with an overall factor A 2 ,/ A 2 = (A r /u c ) . 
Thus it is nonuniversal in the sense discussed above. 

Now return to the correlation function C x (t) which depends on the generalized influence 
functional JF G . As shown above, this can be expressed in the standard form fl4.14p if we 
substitute the modified paths 

£(*') = -©(*'-*)] 
vV) = v(f)-Z[eW)-QV-t)] • 

The effects of the subtractions in Eq. ( |4.21| ) are directly seen in the charge picture. Taking 
into account the constraints ( |4.8|) and (4^), one gets the following changes: In the path 
£(t'), the two blip-charges at times t' = and t' = t are cancelled. In the path fj(t'), the 
sojourn- charge originally located at time t' = is moved to time t' = t. It turns out that 
the influence functionals for the paths of group A and group B are different. We write 



(4.21) 
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n,m 



G H 

^ n,m n,m ' 



*'n,m *- 7 n,m n,m ' \ ) 

A /R 

The blip-interaction-factors G n ' m differ from the standard G n ,m by the absence of the two 
blip charges at t' = and f = t. For group A, this is 

n+m n+m j — 1 



G 



A 

n,m 



exp 



i=i 

j^n+l 



3=2 k=l 



(4.23) 



where Aj k describes the interblip correlations for the modified sequence of charges. If j, k ^ 
n + 1 and ^ n + m, A^ k is again given by ( |4.19| ). In all other cases, the interactions of 
the missing charges have to be dropped in ( [4. 19 ). For instance, for j — n + 1, we have 
^n+i,k = S 2n +2,2k-i - S 2n +2,2k- Similarly, we obtain for group B 

n+m 

G n,m = eX P 



n+m j — 1 
3=2 k=l 



li 

3,k 



(4.24) 
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with analogous modifications in A^ fc for j, k = n and n + m. For instance, we have 
S2n-i,2k — "S , 2n-i,2fe-i- The modified phase factors H l nm take the form 



n.k 



n+m— 1 



H n,m = eX P ™K ^k+l 



k=0 

k^n 

n+m— 1 



(4.25) 



H n,m = eX P ™ K (^k+l + Vn(£n+1 ~ £ 



k=0 

k^n 



Thus each sojourn interacts with the subsequent blip except for sojourn n. For group A, 
the sojourn n is effectively noninteracting, whereas for group B, it effectively interacts both 
with blip n and blip n + m. 

At this point, let us briefly reflect what we have gained so far. First of all, since the 
sequence of the remaining 2n + 2m — 2 blip charges is neutral and comes with a factor 
A 2 ( n+m-1 \ the quantities A and uj c are combined into a factor Ar 2K ^ n+m x \ Thus, the a x 
autocorrelation function turns out to be universal. There is, however, an essential difference 
between the two groups. For group A, the charges in the negative and positive time branches 
are neutral individually. For group B, there is an excess charge ±1 in each branch, and 
only the combined arrangement is neutral again. Since the asymptotic decay of equilibrium 
correlation functions crucially depends on the interactions between the negative and positive 
time branches, we should expect different behaviors for group A and group B. 

Collecting the various results, we obtain explicit expressions for the propagating functions 
of group A and group B in Eq. fl4.lUp , 



Jg(6,*;6>0 + ;77= -£,0 ;r] ,t c 



-■»«££ 

n=0 m=l 



A 2 



2 \ n+m— 1 pt 



/ T^2n,2m-2{tj} ^n.mAi.m ^ 

*° fe} A H} A 



(4.26) 



j g(^^V = 6,0 + ; -£,0 ; 770, to 

00 00 / A2 \ n+m— 1 r-t 



( /\2\ n + m ^ 1 ft 

^EE / ^n-l,2m-l{t,} Yl G lmDn, m £ H% n . (4.27) 

n=lm=n ' Jt ° {^}B { Vj }B 



The summation is over all £j, m = ±1. The superscripts {. . . } A and {. . . } B indicate the 
constraints ( |4.8|) and ( |4.9| ) with £„ +m = respectively. Using Eqs. (f4.1U|) , C x (t) is obtained. 



It is now straightforward to perform the ^-summations and to use symmetry relations under 
exchange {^} — > {— Taking the limit to ~^ —oo, the correlation function becomes inde- 
pendent of the initial value 770- In the end, we find for the symmetrized correlation function 
S x (t) = ReC x (t) and the response function Xx(t) = (— 2/7i) O(t) ImC^t) the expressions 
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1 « -t 

= \ E (-A 2 )" 1 " 1 / 2?o^- a {*i} E <n<2 , (4-28) 

m=l •/-oc {&} A 

OO OO 

= " E E (-A 2 )"^- 1 sin 2 (^) / V 2n . 1>2m ^{ tj } E 6 £n+™ > 

n=l m=2 «/-oo 

£n — £n + 1 — — ^ n + m 

(4.29) 

^(*) = * E E (-A 2 )"^- 1 tan(TrK) / P 2n , 2m _ 2 -fo} E 6 < m C - ( 4 - 30 ) 

n=l m=l {^}A 
„=1 m=l 

x E ^ G lm D £l {sm 2 (7r7r) £ n+1 + cos 2 ^) £ n+m } . (4.31) 
fe} B 

Here we have introduced A 2 = A 2 cos(nK)/2 and 

n+m 

^fti = cos ( e E ^ r i) • ( 4 - 32 ) 



Equations ( |4.28D -( |4.32j ) are exact formal series expansions for the symmetrized equilibrium 
correlation function S x (t) and the response function Xx{t)- Despite their formidable appear- 
ance, we can obtain exact results in certain limits. This is discussed in the remainder of this 
work. For the subsequent analysis, it is convenient to switch from integrations over the flip 
times, Eq. (|4.11|) , to integrations over blip lengths r„- and sojourn lengths s,-. 

V. THE CASE K = \ 

For the value K = |, the above series for S x (t) and Xx(t) can be summed in analytical 
form using the concept of collapsed blips and collapsed sojourns ||. Putting K = | — k with 
k <C 1, the phase factor cos(7rf^) « ttk vanishes in the limit k —>■ 0. In order to have a finite 
contribution for K = ~, each factor cos(tiK) has to be compensated by a 1/k singularity 
arising from the "short-distance" singularity of the breathing mode integral of a dipole (blip 
or sojourn) with interaction e~ s ^ ~ (uj c t)^^ 2k \ Thus we have 

I{K=\)= lim A 2 cos(vrir) / dre" s<T) = - — = 7 . (5.1) 

We shall refer to an expression of the form (|5.1|) as a collapsed dipole. Since a collapsed dipole 



has zero dipole moment, it does not interact with other charges. Further, it is insensitive to 
a symmetric bias factor. In contrast, an odd bias factor in Eq. ( |5.1| ) prevents a dipole from 
collapsing, and combined with a factor cos^-fT), this term vanishes as K — > |. 
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A blip or a sojourn becomes extended when the cos(irK) factor and the short- distance sin- 
gularity are absent. Within an extended blip of length r, the system may make any number 
of visits of duration zero to a sojourn and then returns to the same blip. Mathematically, 
this is described by the insertion of a grand-canonical ensemble of noninteracting collapsed 
sojourns (CS), yielding a CS form factor e~ 7T / 2 . Likewise, within a sojourn of length s, the 
system may make any number of visits of duration zero to a blip state. This is represented 
by a grand-canonical ensemble of noninteracting collapsed blips (CB). Since the system can 
return to either sojourn state, there is a multiplicity factor 2, yielding a CB form factor e~ ls . 

An extended sojourn, say s/-, remains free of insertions only if the subsequent blip is 
weighted with an unconstrained factor in this case, the summation leads to can- 
cellations among the interactions stretching over the extended sojourn, and thus it remains 
bare. It turns out that this is a general rule also for K ^ |, referred to as the £-rule in 
the sequel. For the correlation functions ( |4.29|) - (|4.31|) , e.g., the initial sojourn starting at to 



remains bare due to the factor £1 in the exact formal expressions. There are no other bare 
intervals in the negative-time branch. Thus, the limit ^ —oo is well-behaved. 

Based on these concepts, we now analyze the various contributions. Consider first the 
symmetrized correlation function. Assigning the cos(7rf^) factors to the collapsing dipoles 
as in Eq. (pMQ, there is one cos(7ri\~) factor left in Eq. ( f4.29| ). Thus the contribution S^(t) 
vanishes linearly with k as K — > |. In S£(t), the system dwells in the initial sojourn state r\ 
until time zero. At this time, it hops into the blip state £ = —r\ where it stays until time t, 
resulting in a factor cos(et). The blip of length t is decorated with a CS form factor. Piecing 
the various components together, we find the damped oscillatory behavior 

S x (t) = S£(t) = cos(et) e" 7t/2 . (5.2) 

The contributions to S^(t) are sketched diagrammatically in Fig. 0. Since only collapsed 
sojourns contribute to S x (t) and the short- distance behavior of the pair interaction is inde- 
pendent of temperature, the expression (|5.2|) is valid at any temperature. 



. _q. . 

t t 

FIG. 2. The diagram describing S x (t). The full and dashed lines represent sojourns and blips, 
respectively. The empty box represents the insertion of a CS form factor within the blip interval. 
The bullets mark transitions which are free of bath correlations because of the modified influence 
functional. 

Consider next the response function. The contribution of group A is sketched in Fig. ^. In 
the negative-time branch, the initial sojourn is followed by an extended blip and an extended 
sojourn state. Both of them are equipped with a CS and CB form factor, respectively. At 
time zero, the system hops back into a blip state and stays there until time t. The extended 
blip is again decorated with a CS form factor. In mathematical terms, we have for t > 

l£(t) = - A 2 sin(et) e" 7 * /2 / dr ds sin(er) e~ s{r) e" 7T/2 e" 7S . (5.3) 
* Jo 
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Now, as shown in Ref. ||, the double integral times the factor A 2 is just P^ = (cr z (t — > oo)). 
In the end, we find 

xt(t) = (2/h)P 00 sm(et)e-^ 2 , (5.4) 

2 fl he \ , N 

P 00 = -Imib(- + — -— +i — , 5.5 

where ip{z) is Euler's digamma function. Thus we find again exponential decay, resulting 
from exponential suppression factors due to collapsed blips or sojourns in each interval. 



-□- 



o 



t 



FIG. 3. The diagrams for Xx(t) (top) and Xx(t) (bottom). The full box represents the insertion 
of a CB form factor. The other symbols are analogous to Fig. ^. The upward and downward spikes 
symbolize the remaining charges. 

Now we turn to the contributions of group B. It is immediately clear that the part of 
Xx(t) resulting from the second term in the curly bracket of Eq. (|4.31|) vanishes as k 2 in the 
limit n — > 0, whereas the first term is nonzero in this limit. Here, the system hops from the 
initial sojourn into a blip at a negative time — r and stays there until time zero, where it 
returns to a sojourn state. At time s it hops again into a blip state and dwells in this state 
until time t. Again, each blip interval is decorated with a CS form factor, as discussed above. 
Because of the factor £ n+ i in Eq. (|4.31|) , however, the extended sojourn in the positive time 
branch is free of collapsed blips. The interacting dipole has length r + s and introduces 
correlations between the negative and positive time branches (see Fig. |3|). Thus we have 



a 2 r°° r 



dr I dse-^ t+T - s)/2 e- sir+s) cos[e(t 



t - s)] . (5.6) 

'o Jo 

Introducing the dipole length r + s as a new integration variable, performing the other 
integrations, and combining the resulting expression with Eq. ( |5.4| ), we obtain 

Xx {t) = (2/n)[an{et)F 1 (t) + cos(et)F 2 (t)] (5.7) 



with the functions 



A 2 poo 

Fi (*) = —/ dr e~ 5(T) sin(er) (e^-^ 2 + 
2 7 Jo 

A 2 roc 

F 2 (t) = — J dr e~ 5(r) cos(er) (e^ 1 *^ 



-j(t+r)/2\ 



-r|/2 _ e - 7 (t+r)/2 



)■ 



(5.8) 
(5.9) 
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For asymptotic times t 3> I/7, we find from Eq. (|5.7|) at zero temperature 



X*(t) 



r 



1 



7r/i 7 2 + 4e 2 7^ 



(5.10) 



The algebraic decay law arises from the contribution of group B. Because of the absence 
of collapsed blips in the sojourn interval s, this interval gets effectively very large, s ~ t 
for t ^> I/7. The 1/t law in Eq. ( |5.10| ) is simply the signature of the bare intra-dipole 
interaction, e~ s ® oc 1/t for K — |. The algebraic law at T = is not only of academic 
interest, since it is also valid at low but finite temperatures in the intermediate time regime 
I/7 <C t <C h/3. In the asymptotic limit t 3> K(3 ^> I/7, we find exponential decay, 



Xx{t) 



16 



7 



%ry rf _|_ 4 e 2 ^ 



J_ e -"i*/ 2 



(5.11) 



where v x = 2ir/h(3 is the lowest bosonic Matsubara frequency. 

Since we have calculated the expressions (5T2) and ( |5.7| ) independently, we are now in a 
position to verify whether they are consistent with the fluctuation-dissipation theorem, 



S x (u) = ftcoth (h(3u/2) x " x (u) 



(5.12) 



Taking the Fourier transform of S x (t) and of Xx(t), we find for the spectral function S x (u) 
and the absorptive part of the dynamical coherence susceptibility Xx{^) 



S x (uj) 

x2H 



7 



7 2 /4 + cj 2 + e 2 



( 7 2 /4 + cu 2 + e 2 ) 2 



7 

— tanh 
n 



2k B T 



4e 2 a; 2 ' 
7 2 /4 + uj 2 + e 2 



(5.13) 



( 7 2 /4 + cj 2 + e 2 ) 2 - 4e 2 u; 2 



This confirms that the FDT is satisfied. Finally, the real part x'xi 00 ) °f th e dynamical 
susceptibility reads (we put e = for simplicity) 



nh 7 2 + 4a; 2 



Re 



A-nk B T 



1 



ib[~ + i 
\2 2vrfc D T 



(5.14) 



The linear static susceptibility x x ^ — Xxi 00 ~^ 0) diverges logarithmically as T — > 0. 

One final remark on the case K — 5 is appropriate. The correlation functions can 
also be calculated in a fermionic representation by exploiting the equivalence of the spin- 
boson model for K — ^ with the Toulouse limit of the anisotropic Kondo or resonance level 
model 0. One finds that the a x correlation function in the resonance level model directly 
corresponds to the a x correlation function of the spin-boson model [19|]. The investigation 
of the spin-boson model is convenient when we depart from the particular case K — |. 



VI. THE CASE K 



A. Expansion around K = A 



In the previous section, we have solved the case K 
concept of collapsed blips and collapsed sojourns. Let us now consider the regime K = \ 



I in analytic form by using the 

— K 
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with k < 1 and perform an expansion around the solution of the correlation function for 
K = ~. For finite k, a dipole is actually no longer collapsed. The basic idea now is to 
develop a ^-expansion by systematically taking into account the finite lengths of the blips 
and sojourns. To put up a general computation scheme, it is essential to split the breathing 



mode integral I(K) given in Eq. (|5.1 ) into a contribution h(K) from the short-length interval 
< t < I/7 and a residual contribution Iz(K) from lengths r > I/7. The inverse time 
scale 7 is selfconsistently determined by the short-time part I\(K). We have 

y^i-^^AV.pr^^g 1 " , (6.1) 

where the factor tik is the remnant of the cos(7rf^) phase factor. The frequency 7 is the 
effective inverse time scale of the problem. The short-time part I\(K), representing either a 
collapsed blip or a collapsed sojourn, can be treated exactly in the same manner as described 
in the previous section. That is, all possible arrangements of collapsed dipoles within an ex- 
tended sojourn and blip of length s and r, respectively, add up to a CB form factor e~^ s and 
CS form factor e~^ T//2 . The strategy is then completed by developing a systematic scheme 
to calculate the contributions from the extended blip and sojourn intervals, Tj, Sj > I/7. 
Since the leading contribution hiK) is of order k, it is natural to set up an expansion of 
the correlation function in the number of extended dipoles. Clearly, this is not a systematic 
expansion in k since every extended dipole also contributes higher-order corrections in k. 
However, the strategy will allow to extract the actual long-time behavior of the correlation 
function. To perform the analysis, it is useful to introduce a diagrammatic picture. A 
generic contribution of order n m is obtained by adding m extended dipoles to the respective 
diagram for K — |. Its structure and asymptotic behavior are essentially determined by the 
following rules which are drawn from the exact formal expressions. 

1. ) Insertion of an extended sojourn into a dressed blip interval leads to the diagram of 
Fig. [| (top), whereas insertion of an extended blip into a dressed sojourn interval is dia- 
grammatically represented in Fig. |] (bottom). 

2. ) An extended dipole with insertion of a CB or CS form factor has an effective length of 
order I/7 and therefore cannot produce algebraic decay of the correlation function. 




FIG. 4. Insertion of extended dipoles according to rule 2. 

3.) An extended dipole which is free of CB and CS form factors has a length r >> I/7. 
Therefore, it is sensitive to the unscreened dipole interaction e~ s ^ oc t~ 2K , and its length 
is eventually limited by the overall length t. 

These rules are in correspondence with the above £-rule. We will now apply them to the 
symmetrized correlation function S x (t) and to the response function Xx(t)- 
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B. Response function Xxif) 

We start the discussion of the response function by considering the path contributions of 
group A, Eq. ( |4.30|) . In order K m , the diagram for Xx {P)i Fig- (top), is supplemented by m 
extended dipoles. They can be arbitrarily distributed among the negative and positive time 
branches using rule 1 (cf. Fig. ^). As a result, each interval (except for the first sojourn) is 
dressed by a CB or CS form factor. Thus, using argument 2, Xx(t) decays exponentially. 

Consider next the contribution to X^{t) from the first term in the curly bracket of 
Eq. (|4.31|) , referred to as Xx 1 ^)- I* 1 or der n m , the diagram in Fig. |3| (bottom) for Xx(t) * s 
modified as follows. There are insertions of m extended sojourns which can be arbitrarily 
distributed among the two blip intervals displayed [rule 1]. Each of them is again confined 
to a length of order I/7. Due to the factor £ n+ i in Eq. ( |4.31 ), however, the initial sojourn 



at positive times remains free of insertions. At times t I/7, the length of this interval is 
therefore effectively t. Employing argument 3, we see that the contribution X^(i) decays 
as e~ s ^ oc t~( l ~ 2K \ This law is generally valid, only the prefactor depends on the number 
of extended dipoles considered. 

Starting with order k 2 , there is also a contribution from the second term in the curly 
bracket of (|4.31|) , called Xx 2 (t)- The diagrams are as for Xx 1 ^)^ a P ar t from the crucial 



difference that the first sojourn in the positive time branch is dressed. This is due to the 
absence of the factor in Xx 2 (J)- Thus, according to rule 2, x^ 2 if) decays exponentially. 



C. Symmetrized correlation function S x (t) 

Consider first S^(t), Eq. (|4.28| ), which has dynamics only in the positive time branch. 
Employing rule 1, we have m extended sojourns in order n m , and each of the blip and sojourn 
intervals is dressed. Thus, S^{t) decays exponentially on the time scale I/7. 

As emphasized in the previous section, the leading contribution to S®(t) is of order k. 
This term is found to be 

A 2 r°° /■* 

S*(t) = -7TK— / dr ds e- s ^ e -^ t+T - s)/2 {e^ s - l} cos[e(t - r - s) ) . (6.2) 
2 Jo Jo 

The reason for the subtraction in the curly bracket is the missing of diagrams without any 
insertions in the sojourn interval s. There is always at least one collapsed blip due to the 
constraint in the ^-summation of expression (|4.29|) . Introducing the length r + s as a new 
integration variable, the other integrations can be performed. With the definition 

a 2 roo 

F 3 (t) = — / dre^ sm(er) (e-^/ 2 - e-"'^ 2 ) , (6.3) 
2 7 Jo 

and with 7 replaced by 7 in the expression ( |5.9| ) for i 7 ^), we find 
S*(t) = ttk I [ cos(et)F 2 (t) + sin(et)F 3 (t) ] 

- ( f dT T e- s(r) e-^ +t )/ 2 cos[e(i - r)] (6.4) 
2 ^ Jo 

+ t ^°°dre- s ^e-^ T+t)/2 cos[e(t - r)] ) j . 
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The dominating contribution for t ^> I/7 comes from the first line in Eq. ( |6.4| ), yielding 

The origin of the algebraic decay for k 7^ is the subtraction term in the curly bracket in 
Eq. Q5.2|) . The analysis shows that the subtraction also appears in all higher orders in k. 
Thus the asymptotic behavior £~( 1-2k ) is generally valid for k 7^ 0. For k — > 0, the prefactor 
of the algebraic decay law vanishes and the decay is exponential [cf. Eq. fl5.2|) 1. 



VII. LONG-TIME BEHAVIOR FOR GENERAL K < 1 

From the structure of the various contributions for K = | — k, we can draw conclusions 
for the long-time behavior of the correlation functions for general K < 1. For K substantially 
different from |, the modifications concern the CB and CS form factors inserted in a given 
interval. Instead of collapsed noninteracting blips and sojourns, we now have extended 
interacting dipoles and it is no longer possible to perform the grand-canonical sum in analytic 
form. However, according to the £-rule, the sequences of charges are grouped into clusters 
which are separated by bare sojourns. Because of the alternating sum of the charges within 
a cluster, the length of a cluster is effectively of order 1/A r , where A r is the renormalized 
frequency ( f2~5D - 

With this being the only essential modification, the asymptotic behaviors of the various 
contributions to the correlation function at times t ^> 1/A, r emerge as follows. In group A, 
there is a single neutral cluster surrounding the origin of the time axis. Hence, both S£{t) 
and Xx (t) decay exponentially. In group B, we have a charged cluster in each time branch, 
satisfying overall neutrality. Since in both branches the initial sojourn is free of insertions, 
the two clusters are near the origin and near t, respectively, and they interact with the 
unscreened charge-charge interaction e~ s ^ oc t~ 2K |22| . This interaction directly determines 



the long-time behavior of S^(t) and Xx(t)- The contributions of group B predominate over 
the exponential contributions of group A for t 3> 1/A r . Thus we have asymptotically 

S x (t) oc e" 5(t) oc t~ 2K , K±\ (7.1) 
Xx (t) oc e" 5 « oc t~ 2K . (7.2) 

Thus, the a x autocorrelation function at T = decays with a power law. The power depends 
on the damping strength. Again, the T = decay laws ([7.1Q and ( |7.2| ) hold also at very 
low temperatures in the intermediate time regime 1/A r t <C h/3. In the asymptotic limit 
t 3> h/3 ^> 1/A r , the correlation functions show exponential decay, 

S x (t) oc e~ K ^ , x*(t) « er Kv * , (7.3) 

where the decay rate is K times the lowest bosonic Matsubara frequency v\ = 2ir/hf3. 

Let us now put the decay law ( |7. 1|) in perspective with the generalized Shiba relation for 
the a z correlation function p|-|8|,|r2"f . In the regime t ^> 1/A r , this relation is expressed as 

S z (t) = -2K [hx® /2] 2 \ 2 . (7.4) 
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In a charge representation for S z (t), the 1/t 2 decay law reflects the dipole-dipole interaction 
between a neutral cluster in the negative-time branch and a neutral cluster in the positive 
time branch. The power of the algebraic interaction is independent of the coupling strength 
and is 2 for Ohmic dissipation. 



Our findings are consistent with the fluctuation-dissipation theorem, Eq. ( 5.12|) . Upon 
Fourier transforming Eq. (|7.2|) , we get Xxi. u ^ 0) oc sgn^)!^! 2 ^" 1 for < K < 1. Using 
the FDT relation ( |5.12| ) for T = 0, we obtain S x (u — > 0), and transforming back to time, we 
find consistency with the law ( |7.1| ) for S x (t). As a by-product, we obtain a useful relation 
directly connecting the prefactors of the expressions ( [7.1| ) and (|7.2|), 

S x (t) = (ft/2) cot(nK) , t > 1/A r . (7.5) 



In lowest order in k, this relation is confirmed by the results ( 5 . 1 0|) and ( |6.5|) . The case K — | 



is special, since the prefactor of the 1/t law for S x (t) vanishes according to the relation ( [7.5| ). 
This is in agreement with the result (|6.5| ) obtained from the direct computation of S x (t). 

The asymptotic decay law ( |7.2| ) leads to a different behavior for the linear (zero bias) 
static susceptibility xi = Jo°° dtXxif) for K below and above ~. For K < |, the slow decay 
of Xx{t) implies that the linear static susceptibility diverges algebraically, 

x (0) K jOK-l ag 

T — > 0. This indicates that the system responds to a coherence inducing perturbation oc a x 
in a nonlinear manner. Interestingly, the regime K < | coincides with the coherence regime 
for the population (u z (t)) at zero bias JH]|JT5| ]. For K > 5, the decay of Xxif) is sufficiently 
fast so that the linear static susceptibility is finite at T = 0. This corresponds to the 
incoherent regime for (a z (t)) at zero bias. The transition from nonlinear to linear response 
at K = | thus reflects the intrinsic coherence properties of the system. These properties of 
the static susceptibility have been confirmed numerically in Ref. Jl5[ . 



In conclusion, we have studied within a real-time approach the equilibrium correlation 
function of the polaron-dressed tunneling or coherence operator in the dissipative two-state 
system. This quantity turns out to be universal in the scaling limit. The elimination of the 
bath modes leads to a modified influence functional that can be recast into the standard 
form at the expense of introducing modified system paths. We have obtained the exact 
formal expressions for the coherence correlations for arbitrary damping strength K, and we 
have presented analytic results for the particular case K = | and for the narrow regime 
K = | — k with k < 1. The long-time behavior is found to be oc t~ 2K for general K < 1, 
reflecting the loss of coherence with increasing damping strength. Generally, an algebraic 
decay law reveals that the initial state of the global system is correlated. In a charge picture, 
the particular decay oc t~ 2K expresses the interaction between two clusters with an excess 
charge of opposite sign. Since they are embedded in the vacuum (£-rule), the interaction is 
unscreened. 
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